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Abstract. —Phylogenetic hypotheses about the Cyatheaceae have been based largely on either 
New World or Old World species. In this study we examine molecular data for species from both 
hemispheres. Chloroplast DNA restriction site variation has been compared among 23 species of 
Cyatheaceae from Papuasia and the Pacific and 61 species from the New World. Parsimony anal¬ 
ysis of the variation documented supports the existence of three major clades, and reveals a cor¬ 
relation between phylogenetic relationship and biogeography. Species from the Neotropics appear 
derived with respect to those from the Paleotropics in the Alsophila clade; in this clade species 
from Papua New Guinea and Lord Howe Island are basal to those from the American tropics. A 
similar relationship is tentatively suggested for the Cyatheo clade, where the sole Old World 
representative, (C. howeana ) from Lord Howe Island, is basal to the remainder ol the clade. This 
finding extends the worldwide distribution of the Cyathea clade across the Pacific. No clear cor¬ 
relation between phylogeny and biogeography of the Sphaeropteris clade was observed based on 
current data. These results suggest that additional sampling from other regions in the Old World 
tropics may further elucidate the evolutionary history of this group of plants. 


The Cyatheaceae (the scaly tree ferns) are a family of approximately 500 
species occurring throughout the wet tropics, particularly in cool montane rain 
forests. Regions that are especially rich in species include the Greater Antilles, 
Central America, the northern part of the Andes (including Venezuela, Colom- 


Madaga 


Most 


Whereas 


geographic ranges of species are, for the most part, fairly well known, the 
ranges of genera are not because there is a lack of agreement on generic bound¬ 
aries. This lack of consensus on generic circumscription is illustrated by stud¬ 
ies of Tryon and Tryon (1982), Holttum and Edwards (1983), and Lellinger 
(1987) in which six, one, and four genera were recognized, respectively. Recent 
studies of chloroplast DNA restriction site variation in New World Cyatheaceae 
(Conant et al., 1994; Stein et al., in press) have demonstrated the existence of 
three well-defined evolutionary lineages (the Alsophila clade, the Cyathea clade, 
and the Sphaeropteris clade). The Alsophila clade is represented in all centers 
of tree fern diversity; however, the clade is most diverse in Malesia. The 
Snhaeronteris clade is also most diverse in Malesia. but is absent from Africa 
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and Madagascar and represented in the neotropics by only eight species, the 
S. horrida group (Tryon, 1971). The distribution of the Cyathea clade is con¬ 
troversial. This clade includes the genera Trichipteris, Cyathea, and Cnemi- 
daria, all of which Tryon and Tryon (1982) believed were restricted to the 
American tropics. Holttum (1982), however, argued that the range of this evo¬ 
lutionary lineage extended beyond the American tropics. He proposed that Cy. 
decurrens (Hook.) Copel. ( Alsophila decurrens Hook., of Tryon, 1970) and its 
allies in the islands of the Western Pacific were more closely related to mem¬ 
bers of the Cyathea clade from Central and South America than to members 
of the Alsophila clade as indicated by Tryon (1970). This would suggest that 
the Cyathea clade is not restricted to the Neotropics, but extends westward 
through the islands of the Pacific to Australia and New Ireland. 

Hypotheses about the phylogeny of the Cyatheaceae have been influenced 
by an Old World (Holttum) or a New World (Tryon) perspective. To bring a 
more global point of view to our ongoing investigation of tree fern phylogeny, 
we have obtained new cpDNA restriction site mutation data from the tree ferns 
ot Papua New Guinea and Lord Howe Island. Combining these data with those 
previously obtained (largely from New World taxa) permits a broader exami¬ 
nation of clade membership and the investigation of biogeographic patterns in 
the group. Specifically, our analysis addresses the following questions: 1) Is 
the Cyathea clade restricted to the New World? 2) Is there a close alliance of 
the S. horrida complex with Papuan species? 3) Are the three major clades 
identified earlier maintained as separate lineages? 4) Are Sphaeropteris and 
Cyathea the derived sister lineages as suggested by earlier analyses? 


Materials and Methods 

Total DNA was extracted from fresh, liquid-nitrogen frozen or silica gel-dried 
leaf tissue (Stein, 1993). The study set included 84 species of Cyatheaceae from 
the Caribbean, Central and South America, Papua New Guinea, Lord Howe 
Island, Australia and New Zealand (Table 1). Two outgroup species, Lopho- 
soria quadripinnata (Lophosoriaceae) and Dicksonia gigantea (Dicksoniaceae), 
both from the American Tropics, were also included. DNAs were digested with 
20 restriction enzymes: Apal, Aval, BamHl, Bell, Bgl I, BsfEII, EcoRV, Hindlll, 
Mlul, Nrul, Pstl, Pvull, Sacl, SacII, Safi, Seal, Smal, Sphl, Stul, and Xhol. 
Electrophoresis, Southern blotting, nick translation, and hybridizations were 
carried out as in Stein (1993). 

Detailed restriction site maps were constructed for A. bryophila, A. salvinii, 
Cy. furfuracea, S. cooperi, L. quadripinnata, and D. gigantea (Conant et al., 
1994; Stein et al., in press) using cloned Adiantum capillus-veneris cpDNA 
fragments (Hasebe and Iwatsuki, 1990). Restriction site variation was assessed 
using eight Adiantum probe groups that combined two or three adjacent 
cpDNA fragments (Table 2) in the hybridization experiments. 

Data were analyzed with the computer program PAUP (Phylogenetic Analysis 
Using Parsimony, Version 3.1.1; Swofford, 1992) on a Macintosh Quadra 950. 
Because the number of species was so large, only heuristic searches could be 
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Table 1. Species of Cyatheaceae included in this study. Lophosoria quadripinnata (Lophosoriaceae) 
and Dickson in gigantea (Dicksoniaceae) are outgroup species. Voucher specimens of each species are 
deposited at the Lyndon State College Herbarium (LSC), and a duplicate set at the Harvard University 
Herbaria (OH). Species placed in the genus Nephelea by Tryon (1970) are designated with an asterisk. 
Plants were collected from Costa Rica (CR), the Dominican Republic (DR), Honduras (HO), Lord Howe 
Island (LHI), Papua New Guinea (PNG), and Puerto Rico (PR). LSCG and QBG are the Lyndon State 
College greenhouse and the Quail Botanical Garden, Encinitas, California, respectively. Voucher specimen 
numbers cited are those of Conant. Two vouchers are cited for Sphaeropteris brunei , because the collection 
number on the DNA sample was lost. All species were collected in duplicate so the DNA had to come 
from one of the two specimens cited. 


Species 

1. AI sop hi la ahbottii (Maxon) R. Tryon 

2. /\. amintae Conant 

3. A. brevipinna (Benth.) R. Tryon 

4. A . brooksii (Maxon) R. Tryon 

5. A. bryophila R. Tryon 

6. A . coactilis (Holtt.) R. Tryon 

7. A. cuspidata (Kunze) Conant* 

8. A. erinaceu (H. Karst.) Conant* 

9. A. ferdinandii R. Tryon 

10. A. firma (Baker) Conant* 

11. /V. foersteri (Rosenst.) R. Tryon 

12. A. fulgens (C. Chr.) Conant* 

13. A. geluensis (Rosenst.) R. Tryon 

14. A. gleichenioides (C. Chr.) R. Tryon 

15. A. hooglandii (Holtt.) R. Tryon 

16. A. im ray ana (Hook.) Conant 

17. A. jimeneziana Conant* 

18. A. minor (D.C. Eaton) R. Tryon 

19. A. mueUeri (Baker) R. Tryon 

20. A. nigrolineata (Holtt.) R. Tryon 

21. A. nigropaleata (Holtt.) R. Tryon 

22. A. polystichoides (Christ) Conant* 

23. A. portoricensis (Kuhn) Conant* 

24. A. pseudomuelleri (Holtt.) R. Tryon 

25. A. sale in ii Hook. 

26. A. urban ii (Brause) R. Tryon 

27. A. woodwardioides (Kaulf.) Conant* 

28. Cue midaria amabilis (Morton) R. Tryon 

29. Cn. choricarpa (Maxon) R. Tryon 

30. Cn. consimilis Stolze 

31. Cn. grandiflora (Wi 1 Id.) Proctor 

32. Cn. horrid a (L.) C. Presl 

33. Cn. karsteniana (Klotzsch) R. Tryon 

34. Cn. mutica (Christ) R. Tryon 

35. Cn. spectabilis (Kunze) R. Tryon 

36. Cvathea andina (H. Karst.) Domin 

37. Cv. a r bo re a (L.) Smith 

& 

38. Cv. caracasana (Klotzsch) Domin 

39. Cv. delgadii Sternb. 

40. Cv. divergens Kunze 

41. Cv. J'ulva (Mart. & Gal.) Fee 

42. Cv. furfu race a Baker 


Abbreviation Source Voucher 


ALabb 

DR 

4305 

ALami 

PR 

4318 

ALbre 

LHI 

4667 

ALbro 

PR 

4310 

ALbry 

PR 

4322 

Alcoa 

PNG 

4589 

A Lc u s 

CR 

4427 

A Leri 

CR 

4386 

ALfer 

LHI 

4666 

ALIir 

HO 

4364 

ALfoe 

PNG 

4646 

ALful 

DR 

4287 

ALgel 

PNG 

4661 

ALgle 

PNG 

4584 

ALhoo 

PNG 

4650 

ALimr 

VE 

4466 

ALjim 

DR 

4287 

A Lin in 

DR 

4297 

ALmue 

PNG 

4582 

ALnil 

PNG 

4636 

ALnip 

PNG 

4647 

ALpol 

CR 

4391 

ALpor 

PR 

4321 

ALpse 

PNC. 

4586 

ALsal 

HO 

4365 

ALurb 

DR 

4301 

ALwoo 

DR 

4312 

CNama 

VE 

4495 

CHcho 

CR 

4414 

CNcon 

VE 

4497 

CNgra 

VE 

4487 

CNhor 

PR 

4343 

CNkar 

VE 

4471 

CNmut 

CR 

4385 

CNspe 

VE 

4482 

CYand 

PR 

4334 

CYarb 

PR 

4345 

CYcar 

PR 

4326 

CYdel 

VE 

4455 

CYdiv 

CR 

4384 

CYful 

CR 

4397 

C Y fu r 

PR 

4325 
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Table 1. Continued. 


Species 

Abbreviation 

Source 

Voucher 

43. Cy. holdridgeana Nisman & L.D. Gomez 

CYhol 

CR 

4406 

44. Cv. howeanci Domin 

CYhow 

LHI 

4665 

45. Cy. multiflora Sm. 

CYmul 

VE 

4425 

46. Cy. parvula (Jenman) Domin 

CYpar 

PR 

4332 

47. Cy. speciosa Willd. 

CYspe 

VE 

4475 

48. Cy. suprastrigosa (Christ) Maxon 

CYsup 

CR 

4402 

49. Cy. tenera (Hook.) T. Moore 

CYten 

PR 

4333 

50. Dicksonia gigantea H. Karst. 

DIgig 

CR 

4400 

51. Lophosoria quadripinnata (J.F. Gmel.) C. Chr. 

LOqua 

PR 

4339 

52. Sphaeropteris aeneifolia (Alderw.) R. Tryon 

SPaen 

PNG 

4578 

53. S. angiensis (Gepp) R. Tryon 

SPang 

PNG 

4629 

54. S. atrox (C. Chr.) R. Tryon 

SPatr 

PNG 

4606 

55. S. brunei (Christ) R. Tryon 

SPbru 

CR 

4388/4392 

56. S. concinna (Baker) R. Tryon 

SPcon 

PNG 

4623 

57. S. cooperi (F. Muell.) R. Tryon 

SPcoo 

LSCG 

4304 

58. S. elongata (Hook.) R. Tryon 

SPelo 

CR 

4410 

59. S. glauca (Blume) R. Tryon 

SPgla 

PNG 

4622 

60. S. horrida (Liebm.) R. Tryon 

SPhor 

HO 

4363 

61. S. insignis (D.C. Eaton) R. Tryon 

SPins 

DR 

4396 

62. S. magna (Copel.) R. Tryon 

SPmag 

PNG 

4630 

63. S. medullaris (G. Forst.) Bernh. 

SPmed 

QBG 

4303 

64. S. parianensis Windisch 

SPpar 

VE 

4496 

65. S. procera (Brause) R. Tryon 

SPpro 

PNG 

4645 

66. S. robust a (Watts) R. Tryon 

SProb 

LHI 

4663 

67. S. senilis (Klotzsch) R. Tryon 

SPsen 

VE 

4479 

68. S. tomentosissima (Copel.) R. Tryon 

SPtom 

PNG 

4581 

69. Trichipteris armata (Sw.) R. Tryon 

TRarm 

PR 

4328 

70. T. borinquena (Maxon) R. Tryon 

TRbor 

PR 

4341 

71. T. cordata (Klotzsch) R. Tryon 

TRcor 

VE 

4473 

72. T. decomposita (H. Karst). R. Tryon 

TRdec 

VE 

4440 

73. T. gibbosa (Klotzsch) Barrington 

TRgib 

VE 

4469 

74. T. kalbreyeri (Baker) R. Tryon 

TRkal 

VE 

4448 

75. T. microdonta (Desv.) R. Tryon 

TRmic 

HO 

4357 

76. T. micropityIla (Klotzsch) R. Tryon 

TRmip 

VE 

4468 

77. T. nigripes (C. Chr.) Barrington 

TRnig 

CR 

4417 

78. T. pauciflora (Kuhn) R. Tryon 

TRpau 

VE 

4445 

79. T. procera (Willd.) R. Tryon 

TRpro 

PR 

4338 

80. T. scabriuscula (Maxon) R. Tryon 

TRsca 

HO 

4346 

81. T. schiediana (C. Presl) R. Tryon 

TRsch 

HO 

4368 

82. T. steyermarkii R. Tryon 

TRste 

VE 

4500 

83. T. stipularis (Christ) R. Tryon 

TRsti 

CR 

4395 

84. T. trichiata (Maxon) R. Tryon 

TRtri 

CR 

4422 

85. T. tryonorum (Riba) R. Tryon 

TRtry 

VE 

4465 

86. T. ursina (Maxon) R. Tryon 

TRurs 

CR 

4423 


executed. Searches were carried out with the SIMPLE, CLOSEST and AS IS 
addition sequences with ten trees held at each step during stepwise addition. 
When simple addition sequences were employed, A. abbottii served as the ref¬ 
erence taxon. Starting trees were also constructed in MacClade (Maddison and 
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Table 2. Chloroplast DNA restriction site mutations scored in 84 species of Cyathcaceac and the outgroup 
taxa Lophosoria quadripinnata (Lophosoriaceae) and Dicksonia gigantea (Dicksoniaceac) based on 20 en¬ 
zymes. No phylogenetically significant mutations were found for the enzymes BstE II, Min I, and PvuU. 
Only synapomorphies are listed. Mutations are numbered according to the alphabetical listing by enzyme 
(No.) and by the actual column in the data matrix (Col. Matrix). Mutations are listed with ancestral state 
lirst and derived state last. Species with mutated DNAs are listed as assigned by PAUL in the tree in Fig. 
1. Mutations were polarized using Lophosoria quadripinnata and Dicksonia gigantea. Probe groups corre¬ 
spond to Pst\ fragments of cpDNA from Adiantum capillis-veneris (Hasebe and Iwatsuki, 1990). 


Col. 


No. 

Matrix 

Enzyme 

Probe Group 


Mutation 

Species with Mutated DNAs 

1. 

58 

A pa I 

A9 + 5 



7.0 + 

10.0 - 17.0 

28-35, 47 

2. 

59 

A pa I 

A10 + 4 

4 - 

17 

13.3 = 

= 5.5 + 7.8 

28-43, 45-49. 58. 64. 67, 









69-86 

3. 

63 

Apa\ 

A19 T 6 


12 + 3 

9.6 = 

4.2 + 5.4 

28-31, 33-35, 39-41, 43, 45, 




+ 20 4 - 

16 



47^48, 58. 64, 67, 71-78. 









80-86 

4. 

56 

A va I 

All 4- 1 



6.9 = 

0.9 + 6.0 

28-^49, 52-86 

5. 

57 

A i j a 1 

A3 + 20 

4 - 

16 + 

2.9 + 

6.4 - 9.3 

28-49. 58. 64. 67. 69-86 




11 + 1 






6. 

41 

BamH\ 

All + 1 



6.5 = 

1.7 + 4.8 

12. 17. 27 

7. 

42 

Bam HI 

All + 1 



6.5 - 

2.3 + 4.2 

28-36, 38-49, 64, 67. 70-79. 









81-82, 86 

8. 

44 

BamH\ 

A9 + 5 



3.9 + 

6.9 - 10.8 

28-29, 32, 34-35, 47 

9. 

43 

Bam HI 

A9 + 5 



9.7 = 

2.8 + 6.9 

9, 13, 15, 21, 28-49. 52-80, 









82-86 

10. 

79 

Bam HI 

A19 + 6 

4 - 

12 

16.3 = 

- 6.3 + 10.0 

73, 82 

11. 

80 

Bam HI 

A19 + 6 

4 - 

12 

16.0 = 

= 3.2 + 12.8 

1-2, 4-5, 7-8, 10. 12. 16-18. 









22-23, 26-27 

12. 

45 

Bam HI 

A19 + 6 

4- 

12 

3.7 - 

0.7 + 3.0 

45, 77, 86 

13. 

81 

BamH\ 

A10 + 4 

4- 

17 

1.35 + 1.35 = 2.7 

58, 80, 83-85 

14. 

51 

Bcl\ 

A2 + 15 



2.5 = 

0.9 + 1.6 

52-57, 59-63, 65-66, 68 

15. 

55 

Bell 

A1 1 + 1 



5.0 + 

9.3 = 14.3 

52-57, 59-63, 65-66, 68-69, 









80, 83-85 

16. 

52 

Bell 

All 4- 1 



9.3 - 

1.0 + 8.3 

1-27 

17. 

53 

Bell 

A9 4- 5 



1.8 + 

4.6 — 6.4 

28-35, 47 

18. 

54 

Bell 

A8 + 13 

4 - 

2 + 15 

11.4 = 

- 3.6 + 7.8 

1-24, 26-27, 44 

19. 

26 

Bgl I 

A2 + 15 



25.0 = 

- 10.0 + 15.0 

69, 80. 83-85 

20. 

25 

Bgll 

A2 + 15 



26.0 = 

= 8.0 + 18.0 

1-2, 4-5. 7-8, 10, 12, 16-18, 









22-23, 26-27 

21. 

85 

Bgll 

A2 + 15 



45.4 = 

- 19.4 + 26.0 

1-2, 4-24, 26-27 

22. 

86 

Bgll 

A8 + 13 



19.6 = 

- 4.8 + 14.8 

13. 15, 21 

23. 

27 

Bgll 

All 4- 1 



28.5 = 

= 13.5 + 15.0 

28-49, 52-86 

24. 

46 

EcoRV 

A19 + 6 

4 - 

12 

7.1 = 

1.7 + 5.4 

1-24, 26-27 

25. 

6 

EcoRV 

A19 4 - 6 

4 - 

12 

9.3 = 

2.2 + 7.1 

1-35. 37-49. 52-86 

26. 

47 

EcoRV 

A8 + 13 



1 1.7 = 

= 5.2 + 6.5 

44, 52-57, 59-63, 65-66, 68 

27. 

48 

EcoRV 

A4 + 10 

4- 

17 

17.7 - 

- 7.2 + 10.5 

58. 69. 80. 83-85 

28. 

49 

Hindi 11 

A19 + 6 

4- 

12 

5.5 = 

1.9 + 3.6 

1-2, 4-5, 12, 17-18. 23, 26-27 

29. 

39 

Hindi 11 

A19 4- 6 

4 - 

12 

1.7 + 

1.8 - 3.5 

52-57. 59-62, 65-66, 68. 79 

30. 

91 

Hindlll 

A8 + 7 + 

14 

13.4 = 

= 4.2 + 9.2 

28-49, 52-86 

31. 

92 

H indlU 

A8 + 7 + 

14 

1.2 + 

9.2 = 10.4 

28-43, 45-49, 53. 58. 63-64. 









67. 69-78. 80-86 

32. 

93 

Hindi 11 

A8 + 7 + 1 

14 + 9 + 

1.2 + 

4.2 = 5.4 

46. 56. 60-61, 79 


5 
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Table 2. Continued. 


Col. 


No. 

Matrix 

Enzyme 

Probe Group 



Mutation 

Species with Mutated DNAs 

33. 

50 

Hindm 

A9 + 5 



0.8 + 

1.0 = 1.8 

1-24, 26-27 

34. 

72 

Nrul 

A8 + 13 



27.0 = 

= 5.5 + 21.5 

1-2, 4-5, 12, 17-18, 23, 









26-27, 51 

35. 

30 

Nru I 

A8 + 13 



26.5 = 

= 10.0 + 16.5 

45, 64. 67, 70-73, 75-77, 79, 









81-82, 86 

36. 

31 

Nrul 

A8 + 13 

+ 2 + 

15 

6.5 4 

27.0 - 33.5 

52-57, 59-63, 65-66. 68 

37. 

32 

Nrul 

A3 + 20 

+ 16 


9.5 - 

2.0 + 7.5 

3, 9. 13-15, 19, 21, 24-25, 
28-49, 52-59, 62-86 

38. 

33 

Nrul 

All 4- 1 



18.2 = 

= 4.2 + 14.0 

52 57, 59 63, 65-66. 68 

39. 

34 

Nrul 

A9 + 5 



0.9 + 

1 1.6 = 12.5 

38-43, 48-49, 74, 78 

40. 

35 

Nrul 

A9 + 5 



15.3 = 

= 1.8 + 13.5 

46, 69, 83-85 

41. 

89 

Nrul 

A9 + 5 



1.0 + 

3.4 - 4.4 

44, 52-54, 57, 59, 62, 66, 68 

42. 

74 

Nrul 

A18 + 7 

+ 14 


15.3 = 

= 5.2 + 10.1 

1 35, 37-44, 47^49, 63, 74. 78 

43. 

1 

Pstl 

A18 + 7 

+ 14 


1.8 + 

6.3 = 8.1 

1-24, 26-27, 51 

44. 

64 

Pstl 

A18 + 7 

4 14 


19.6 = 

= 6.6 + 13.0 

28, 30-33, 35, 40-41, 47 

45. 

65 

Pstl 

A18 + 7 

+ 14 


19.6 = 

= 8.1 + 11.5 

1-2, 4-5, 12, 17-18, 23, 26-27 

46. 

66 

Pstl 

A18 + 7 

4 14 


0.6 4 

21.4 = 22.0 

52-57, 59-63, 65-66, 68 

47. 

3 

Pstl 

A3 + 20 

4 16 


7.9 + 

8.1 = 16.0 

1-2, 4-5, 7-8, 10-12, 16-18, 









22-23, 26-27, 37, 55 

48. 

67 

Pstl 

A3 4- 20 

+ 16 


6.8 + 

7.4 = 14.2 

1-27 

49. 

2 

Pstl 

A9 + 5 



8.0 = 

0.6 + 7.4 

28, 30-49, 58. 60, 64, 67, 









69-82, 84-86 

50. 

19 

Sad 

A3 + 20 

+ 16 


22.0 = 

= 9.0 + 13.0 

1-2, 4-8, 10-12, 14, 16-20, 









22-24, 26-27 

51. 

20 

Sacl 

All 4- 1 



13.7 = 

= 6.7 + 7.0 

28—43, 45-49, 58, 64, 67, 69-86 

52. 

22 

Sad 

A9 + 5 



8.8 = 

1.3 + 7.5 

52-57, 59-63, 65-66, 68 

53. 

21 

Sad 

A9 + 5 



5.7 + 

8.7 = 14.4 

1-49, 52-86 

54. 

23 

Sad 

A2 + 15 



3.4 4 

24.0 = 27.4 

1-24, 26-27, 64, 67, 71-72 

55. 

24 

Sad 

A8 + 13 



26.0 = 

= 1 1.0 + 15.0 

29, 34 

56. 

87 

SacII 

A9 + 5 + 10 + 

4 

13.0 = 

= 6.4 + 6.6 

9, 13, 15, 21 




+ 17 






57. 

28 

Sadi 

A10 + 4 

+ 17 


12.0 + 17.0 = 29.0 

69, 80, 83-85 

58. 

29 

Sacl I 

A1 1 + 1 



14.1 = 

= 6.5 + 7.6 

1-24, 26-27 

59. 

73 

Sacl I 

A1 1 + 1 

4 18 4 7 

3.2 + 

14.1 = 17.3 

37, 69, 80, 83-85 




+ 14 






60. 

10 

Sail 

A10 + 4 

+ 17 


33.0 = 

= 13.0 + 20.0 

29, 32, 34-37, 42-43, 45, 47. 









58, 64. 67, 69-73, 75-86 

61. 

84 

Sail 

A19 + 6 

+ 12 


4.6 + 

10.5 = 15.1 

31. 52-57, 59-63, 65-66, 68 

62. 

11 

Sail 

All H- 1 



19.8 = 

= 5.8 + 14.0 

1-15, 17-24, 26-49, 58, 64, 









67, 69-86 

63. 

9 

Sail 

A9 + 5 



4.8 = 

2.4 + 2.4 

28-49, 58, 64, 67, 69-86 

64. 

36 

Sea I 

A8 + 13 



10.4 = 

= 3.0 + 7.4 

1-2, 4-5, 7-8, 10, 12-13, 









15 18,21 23.26-27 

65. 

90 

Sea I 

A8 + 13 



18.6 = 

= 8.0 + 10.6 

60-61 

66. 

77 

Sea I 

A8 + 13 

4 2 4 

15 

18.6 = 

= 3.2 + 15.4 

1-2, 4-5, 12, 17-18, 23, 26-27 

67. 

78 

Sea I 

A8 + 13 

4 2 4 

15 

18.6 = 

= 5.5 + 13.1 

28-43, 45-49, 58, 64, 67, 









69-86 

68. 

75 

Sea I 

A8 + 13 

+ 2 + 

15 

1.0 + 

10.4 = 1 1.4 

38-41, 43, 48-49, 74, 78 

69. 

76 

Sea I 

A8 + 13 

4 2 4 

15 

13.1 = 

= 1.1 + 12.0 

28-37, 42^43, 45, 47-48. 58, 


64, 67, 69, 71-73, 75-86 
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Table 2. Continued. 


Col. 


No. 

Matrix 

Enzyme 

Probe Group 


Mutation 

Species with Mutated DNAs 

70. 

37 

Sea 1 

A19 + 6 

+ 

12 

14.2 = 

= 3.6 + 10.6 

28-49, 58, 64, 67, 69-86 

71. 

38 

Sca\ 

Al 1 + 1 



25.0 = 

= 8.0 + 17.0 

1-27, 52-57, 59-63, 65-66, 68 

72. 

40 

Sea 1 

A9 + 5 



3.1 - 

1.3 + 1.8 

1-2, 4-5, 12, 16-18, 23, 26-27 

73. 

61 

Sma I 

Al 1 + 1 



22.5 = 

= 3.5 + 19.0 

28 31, 33 35, 37, 44-46, 64, 67, 









70-73, 75-77, 79. 81-82, 86 

74. 

62 

Smal 

Al 1 + 1 



21.8 = 

- 3.8 + 18.0 

45, 58, 77, 86 

75. 

60 

Smal 

Al 1 + 1 



22 5 = 

= 1.5 + 21.0 

58. 69, 80, 83-85 

76. 

18 

Sph I 

Al 1 + 1 



11.5 = 

= 3.0 + 8.5 

28-49, 58, 64, 67, 69-86 

77. 

12 

S t u I 

A8 + 13 



23.0 = 

= 9.0 + 14.0 

9, 45, 77, 86 

78. 

13 

Stu I 

A10 + 4 

+ 

17 

8.5 = 

2.1 + 6.4 

28-49, 58, 64, 67, 69-86 

79. 

14 

Stul 

Al 1 + 1 



19.0 = 

= 6.3 + 12.7 

38-41, 49, 74 

80. 

69 

Slid 

Al 1 + 1 



18.6 = 

= 5.6 + 13.0 

1-2, 4-5. 12, 17-18. 23, 









26-27, 45, 64, 67. 71-73, 
75-77, 79, 81-82, 86 

81. 

70 

Slid 

Al 1 + 1 



19.0 = 

= 2.4 + 16.6 

28-35, 46-47 

82. 

15 

Sr id 

A9 + 5 



2.1 + 

5.0 = 7.1 

45-46, 77. 86 

83. 

16 

Slid 

A3 + 20 

+ 

16 

1.8 + 

4.4 = 6.2 

28-49, 58 , 64, 67 , 69-86 

84. 

71 

Stul 

A3 + 20 

4" 

16 

0.6 + 

7.8 = 8.4 

37. 58, 69, 80. 83-85 

85. 

17 

Slid 

A19 + 6 

+ 

12 

17.6 = 

= 2.6 + 15.0 

1-2, 4-5, 7-8, 10, 12, 16-18, 









22-23, 26-27 

86. 

4 

Xhol 

A8 + 13 



11.3 = 

= 3.8 + 7.5 

28-49, 52-86 

87. 

5 

Xhol 

A3 + 20 

4" 

16 

8.8 - 

1.0 + 7.8 

29. 32, 34. 36, 43, 45, 47, 57, 
64, 67, 70-73, 75-77, 79, 
81-82, 86 

88. 

7 

Xho\ 

A19 + 6 

+ 

12 

16.9 = 

= 8.3 + 8.6 

25, 79 

89. 

8 

Xhol 

A 19 + 6 

+ 

12 

16.8 = 

= 6.8 + 10.0 

58, 69, 80, 83-85 

90. 

68 

Xhol 

A9 + 5 



0.7 + 

6.0 = 6.7 

29. 69, 80. 83-85 

91. 

82 

Xhol 

A9 + 5 



10.3 = 

= 3.5 + 6.8 

9, 13, 15, 21 

92. 

83 

Xhol 

A 9 + 5 



2.1 4- 

10.3 = 12.4 

1-2, 4-5, 7-8, 10. 12, 16-18, 


22-23, 26-49, 58, 64, 67, 
69-86 


Maddison, 1992) and then imported into PAUP for branch swapping. The Tree 
Bisection Reconnection (TBR) algorithm was used for branch swapping on each 

of the starting trees. The MULPARS, COLLAPSE ZERO-LENGTH BRANCHES, 
and STEEPEST DESCENT options were in effect with MAXTREES set at 20,000. 
Majority-rule consensus trees were computed for each set of 20,000 trees and 
the trees compared for common monophyletic groups. All character reconstruc¬ 
tions were obtained using the ACCTRAN character-state optimization option. 
All searches were carried out using Wagner parsimony with the constraint that 
the ingroup is a monophyletic group (Stein et al., in press). 

Random addition sequences were also used to obtain a broader distribution 
of starting trees upon which to carry out branch swapping with TBR (Maddison 
et ah, 1992). Up to five equally parsimonious trees were saved at the end of 
each replication. After 100 replications, majority-rule consensus trees were com¬ 
puted from all trees stored in memory. Twenty-five searches (100 replications 
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each) were carried out. The frequencies of various monophyletic groups in each 
search were tabulated, and mean frequencies of these groups computed. 

The degree of support for monophyletic groups was also evaluated via boot¬ 
strapping (Felsenstein, 1985; Sandersen, 1989). Again, because of the number 
of taxa, only heuristic searches could be carried out. We used the AS IS ad¬ 
dition sequence, and held one tree at each step during stepwise addition. The 
TBR algorithm was used to carry out branch swapping. The STEEPEST DE¬ 
SCENT and MULPARS options were not in effect. At the end of 100 replica¬ 
tions, a majority rule consensus tree was constructed. 

Tryon and Tryon (1982), Nishida (1989), and Conant et al. (1994) had indi¬ 
cated that Lophosoria and Dicksonia are the closest relatives of the Cyathea- 
ceae, and noted that it is uncertain which of these taxa is closer to the scaly 
tree ferns. Therefore, we have used both taxa as outgroups in polarizing the 
mutations in our data set and compared searches using both species to those 
in which either was used as the outgroup. 

The lengths of trees with alternative arrangements of the major clades were 
detfirmirmri usinv the constraints ontion in PAUP and with MacClade. 


Results 

In this study we have compared cpDNA restriction site variation in 84 spe¬ 
cies of Cyatheaceae. Sixty-one of these taxa are from the American tropics, 
representing about 33% of the species diversity in the New World (Tryon and 
Tryon, 1982). Twenty-three species were from Malesia and the Islands of the 
Western Pacific, accounting for about 10% of the diversity from that region 
(Holttum, 1963; 1964). From this comparison, we found 92 variable restriction 
sites that were phylogenetically informative (Table 2). 

The shortest trees found were 160 steps long. Because the shortest tree ob¬ 
tained is dependent on the starting tree, we used a variety of different heuristic 
search procedures in analyzing the data. This included searches where 20,000 
trees were generated from a single starting tree and searches that involved 
many starting trees derived from random addition sequences. Random addi¬ 
tion sequences enabled us to obtain a wide variety of starting trees (2500) in 
order to broadly sample tree space. Overall, more than 32,767 trees (the max¬ 
imum number that can be stored in memory in PAUP) of 160 steps were found. 
This is not unexpected, because there were only 92 characters available to 
resolve relationships among 84 taxa. When a monophyletic group was present 
in 95% or more of all minimal trees observed, we considered it a feature of 
the set of most parsimonious trees that merited discussion. These monophy¬ 
letic groups are supported by one or more synapomorphies. 

We have also carried out bootstrap analyses of our data set in order to obtain 
another measure of the degree of support for monophyletic groups. In the dis¬ 
cussions below we have referred to the mean frequency among consensus trees, 
the number of synapomorphies supporting a group, and the bootstrapping val¬ 
ues for various monophyletic lineages. All of these parameters are considered 
together to gain an understanding of the degree of support for a group. 
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With the addition of the Old World taxa, the data continue to support the 
existence of three major evolutionary lineages in the family: The Alsophila , 
Cyatheci, and Sphaeropteris clades (Fig. 1). These lineages are present in 100% 
of the consensus trees and are supported by 7, 12, and 7 synapomorphies, 
with bootstrapping values of 80, 62, and 85%, respectively. The results re¬ 
ported here are based on searches in which both Lophosoria and Dicksonia 
were used as outgroups (Fig. 1), but were very similar to those using either 
taxon to polarize the mutations. 


Within the Alsophila clade, the most derived group consists of ten species 
from the Caribbean. This group is present in 100% of the consensus trees, is 
supported by 5 synapomorphies, and has a bootstrapping value of 100%. 
These species, together with other American Alsophila (with the exception of 
A. salvinii), form a larger monophyletic group that is present in 100% of the 
consensus trees, is supported by 5 synapomorphies, and has a bootstrapping 
value of 99%. American Alsophila is derived with respect to other Alsophila 


species analyzed from Papua New Guinea and Lord Howe Island. Together, 
these American and Old World species form a group which is present in 100% 
of the consensus trees, supported by 7 synapomorphies, and occurs in 80% of 
bootstrap replicates. However, A. salvinii, from southern Mexico and Central 
America, lies outside this clade. Alsophila salvinii is linked to the base of the 
Alsophila clade by two mutations. We are currently making plans to collect 
putative relatives of A. salvinii in Borneo to clarify the relationship of this 
species to other members of the family. 

Tryon and Tryon (1982) argued that the taxa within the Cyathea clade were 
restricted to the American tropics and that Gy. decurrens and its allies in the 
western Pacific were members of the genus Alsophila. Holttum (1982), how¬ 
ever, suggested that the latter species were related to Cy. multiflora and its 
allies in Central and South America. In our analysis, Cy. howeana, a member 
of the Cy. decurrens group endemic to Lord Howe Island, clusters with the 
C^yathea clade in 100% of the consensus trees, shares 12 synapomorphies with 
other members of this clade, but is present in only 62% of the bootstrap rep¬ 
licates. If Cy. howeana is indeed a member of the Cyathea clade as hypothe¬ 


sized by Holttum (1982), rather than a member of the Alsophila clade as sug¬ 
gested by Tryon and Tryon (1982), this would extend the geographic distri¬ 
bution ol the Cyathea clade through the islands of the Pacific. However, mov¬ 
ing Cy. howeana to the Alsophila clade increases tree length by only one step. 
Thus the data only weakly suggest that Cy. howeana is a member of the Cy¬ 


athea clade. Additional data might better resolve this problem, but it is inter¬ 
esting to note that in the shortest trees Cy. howeana is basal to all New World 
members of the Cyathea clade. Thus it might be that, as in the Alsophila clade, 
the New World species are more derived than those of the Old World. 

In the Sphaeropteris clade, there are 14 species belonging to Sphaeropteris 
subgenus Sphaeropteris section Sphaeropteris. Within this section, the cpDNA 
data provide evidence for two monophyletic groups. One, including two 
American species, S. horrida and S. insignis, plus S. concinna from Papua 
New Guinea, is present in 94% of the consensus trees, is supported by 1 syn- 
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apomorphy, and has a bootstrapping value of 62%. Tryon (1971), in a study 
of the American tree ferns allied to S. horrida, stated that, “The relation of the 
S. horrida group to species such as S. medullaris of New Zealand and S. con- 
cinna of New Guinea is so close that it must be accepted as a true (evolution¬ 
ary) neotropic-paleotropic relationship.” The cpDNA data support a close re¬ 
lationship of S. horrida and S. insignis to S. concinna, confirming, at least in 
part, Tryon’s prediction. Sphaeropteris brunei, a member of the S. horrida 
group from Central America, appears to lie outside this group, as it lacks one 
mutation present in S. horrida, S. insignis, and S. concinna. However, S. med¬ 
ullaris, the most basal member of the Sphaeropteris clade, does not appear to 
be a close relative of the S. horrida group, as suggested by Tryon. The other 
monophyletic group, including 6 species from Papua New Guinea, S. cooperi 


Queensland 


synapomor 


ping value of 52%. Although members of this group share only one synapo- 
morphy, the cohesive geographical distribution suggests that it may be a nat¬ 
ural evolutionary lineage restricted to Papuasia. 

It is important to determine the relationships of the major clades to one 
another in order to resolve long standing arguments about the direction of 
morphological character evolution in the family. Conant et al. (1994), in a 
study of 23 species of Cyatheaceae, reported that trees in which the Cyathea 
and Sphaeropteris clades were positioned as derived sister groups were four 
to five steps shorter than alternative arrangements of the clades. Stein et al. 
(in press), in an analysis of 63 tree fern species, also found that trees with the 
Cyathea and Sphaeropteris clades as sister groups were more parsimonious 
than other arrangements, although in their study, the difference was only one 


matrix has undermined 


major 


arrangement 


We 


acters that will help to resolve the relationships of the three major clades. 

In summary, we find that in two of the three major clades of the Cyatheaceae 
the basal-most members are from the Old World or are New World species that 
belong to a primarily Old World group (such as A. salvinii). It is premature to 
reach conclusions about the center of origin of the family because we have not 
yet included species from several floristically important regions of the paleo- 
tropics. The evidence available from the fossil record, particularly from the 
study of spores (Mohr and Lazarus, in press), suggests that the derived genus 
Cnemidaria, now restricted to the Neotropics, once had a much broader distri¬ 
bution in Antarctica, Australia, and Asia. Thus, current geographical distribu- 


taxa 


lm 


history of the Cyatheaceae depends on an appreciation of their phylogeny. Sig¬ 
nificant insights into the evolution of this group will only come from data sets 
representing both hemispheres, such as the results reported here. 
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Alsophila Clade 




American Alsophila Old World Alsophila Neotropical 
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Fig. 1. A representative tree derived from an analysis of chloroplast DNA restriction site muta¬ 
tions in the Cyatheaceae. The length is 160 steps (excluding autapomorphies). The consistency 
index is 0.575, and the retention index is 0.939. The major clades and monophyletic groups dis¬ 
cussed in the text are indicated with brackets at right. Species with asterisks in the Alsophila 
clade were included in the genus Nephelea by Tryon and Tryon (1982). The numbers above the 
nodes indicate the number of restriction site mutations. Numbers below the nodes are the mean 
frequency of monophyletic groups in 25 heuristic searches (each of 100 replications) using random 
addition sequences and bootstrapping values (in parentheses). See Table 1 for list of species names 
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that correspond to abbreviations. The tree was constrained so that all Cyatheaceae were together 
and outgroup species were excluded from the ingroup. 
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